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Magnetic-Field-Induced insulator-conductor transition 
in SU (2) quenched lattice gauge theory 
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We study the correlator of two vector currents in quenched SU (2) lattice gauge theory with a 
chirally invariant lattice Dirac operator with a constant external magnetic field. It is found that 
in the confinement phase the correlator of the components of the current parallel to the magnetic 
field decays much slower than in the absence of a magnetic field, while for other components the 
correlation length slightly decreases. We apply the maximal entropy method to extract the corre- 
sponding spectral function. In the limit of zero frequency this spectral function yields the electric 
conductivity of the quenched theory. We find that in the confinement phase the external magnetic 
field induces nonzero electric conductivity along the direction of the field, transforming the system 
from an insulator into an anisotropic conductor. In the deconfinement phase the conductivity does 
not exhibit any sizable dependence on the magnetic field. 

PACS numbers: 11.30.Rd; 12.38.Gc; 13.40.-f 



Recently, heavy ion experiments at the BNL Relativis- 
tic Heavy Ion Collider (RHIC) have found an evidence 
[1] for the so-called chiral magnetic effect [2] in quark- 
gluon plasma. The essence of the effect is the genera- 
tion of electric current along the direction of the external 
magnetic field in the background of topologically non- 
trivial gauge field configurations. Experimentally, the 
effect manifests itself as the dynamical enhancement of 
fluctuations in the numbers of charged hadrons emitted 
above and below the reaction plane in off-central heavy- 
ion collisions. Recently this effect has been studied also 
in lattice gauge theory, and the evidence for charge sepa- 
ration in magnetic field has been found [3, 5]. In [3] it has 
been found that the fluctuations of the electric current 
along the magnetic field are strongly enhanced as com- 
pared to the fluctuations of current in the perpendicular 
directions. This conclusion was also confirmed by an an- 
alytical calculation in the instanton gas model [6]. The 
result of [3] on the difference of longitudinal and trans- 
verse electric current susceptibilities has been reproduced 
later by an analytical calculation [7]; the frequency de- 
pendence of the conductivity has also been evaluated - 
for the weak coupling result, see [8]. 

A natural question to ask is whether this enhance- 
ment of current fluctuations corresponds to a real flow 
of charge, or is just caused by short-lived quantum fluc- 
tuations. This question can be answered by studying 
the current-current correlation functions. The currents 
which correspond to a real transport of charged parti- 
cles should have long-range correlations in time, while 
quantum fluctuations are typically characterized by a fi- 



nite correlation time [9]. Recalling Green-Kubo relations, 
one can see that this property is intimately related to 
the electric conductivity - namely, the real transport of 
charged particles can occur only in conducting media. In 
this Letter we study the tensor of electric conductivity of 
the vacuum of quenched SU (2) lattice gauge theory in 
external magnetic field. We find that the magnetic field 
induces nonzero electric conductivity along its direction, 
transforming the confining vacuum from an insulator into 
an anisotropic conductor. 

Electric conductivity can be extracted from the corre- 
lator of two vector currents ji (x) = q (x) "fiq (x): 



GV,(t)= / cPx{jS^)ij(x,T)) 



(1) 



Following [10], let us define the spectral function p(w) 
which corresponds to the correlator (1) 



Gij (t) = —K (w, t) p i:j (i 
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where T is the temperature. The Kubo formula for the 
electric conductivity then reads [9, 10]: 
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In the limit of the weak time-independent electric field 
Ek, one has (ji) = (JikEk- Thus electric conductivity is 
related to the behavior of the spectral function at small 



2 



frequencies. If there is a gap in the spectrum so that 
Pij (w) — for w < w c , electric conductivity is zero and 
Gij (t) ~ cosh (ui c (t — 2^)) ■ On the other hand, if p (w) 
is not zero near w = 0, one can expect slow nonexponen- 
tial decay of (r). 

To measure the correlator (1), we perform lattice 
Monte-Carlo simulations of quenched SU (2) lattice 
gauge theory. Since quark chirality is very important 
for magnetic effects in non-Abclian gauge theories [2], 
we use the overlap lattice Dirac operator V with exact 
chiral symmetry [11] to measure the vector currents. We 
consider the two-current correlator in the meson chan- 
nel, which is represented in terms of Dirac propagators 
in fixed Abclian and non-Abclian gauge fields and is then 
averaged over an equilibrium ensemble of non-Abelian 
gauge fields A^: 



q (x) 7 ?; g (x) q (y) jjq (y) ) 

Sym [A/j] ^j- 
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where Sym [A^] is the lattice action for gluons A^. A 
uniform magnetic field is added to the Dirac operator 
by substituting su (2)-valued vector potential A^ with 
u (2)-valued one A^ij — > A^ij + 1/2 F^ v x„5ij. In order 
to account for periodic boundary conditions we introduce 
an additional twist for fermions [3, 4]. The quark mass is 
fixed in lattice units at a small value am = 0.01. Previous 
studies of mesonic correlation functions with an overlap 
Dirac operator indicate that the vector current correlator 
depends very weakly on quark mass [12]. 

Strictly speaking, the correlator (5) corresponds to the 
correlator of charged currents, for example u^^d. The 
correlator of neutral currents considered in our Letter, 
j M = c?7 M d, should also contain the disconnected part. 
This part is quite intricate for an accurate numerical 
treatment. We have roughly estimated its contribution 
by inverting the Dirac operator on a subspace spanned 
on some small number M ~ 30 of the lowest Dirac eigen- 
modes, as in [1.7]. It turned out that this part of the 
full neutral current correlator behaves similarly to the 
connected one (5). We do not reproduce these estimates 
here due to uncontrollable systematic errors [17]. 

We use the tadpole-improved Wilson-Symanzik action 
(see, e.g., Eq. (1) in [13]). For inversion, we use a Gaus- 
sian source with radius r = 1.0 in lattice units in both 
spatial and time directions and a point sink (that is, 
quark position is smeared over a Gaussian profile). Wc 
have found that such smearing significantly improves the 
convergence of the maximal entropy method [10, 15, 16] 
at small lattice sizes, while the value of the conductiv- 
ity is practically unaffected. Our lattice parameters are 
summarized in Table I. A uniform magnetic field is in- 
troduced into the Dirac operator as described in [3]. In 
order to obtain the Dirac propagator, wc implement the 
shifted unitary minimal residue method of Ref. [18]. 
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TABLE I. Lattice parameters used in our simulations. The 
critical temperature of the deconfinement phase transition in 
quenched SU (2) gauge theory is T c = 313.(3) MeV [14]. 



It is clear that since the magnetic field is parallel to the 
z axis, the principal axes of the tensor a irj (r) will be the 
x, y and z axes and it is sufficient to consider only the 
diagonal components an (no summation over i=x, y, z). 

We plot some correlators at different temperatures and 
magnetic fields on Fig. 1. The data are for the 14 4 lat- 
tice with spacing a = 0.102 fm (left) and for the 16 3 x 6 
lattice with spacing a = 0.095 fm (right). For the lat- 
ter lattice the temperature is T = 350 MeV = 1.12 T c 
and the theory is in the deconfinement phase. In the 
quenched theory the critical temperature of the decon- 
finement transition is not affected by the magnetic field. 
The temperature T = 1.12 corresponds to the chirally 
restored phase. 

One can see that without the magnetic field the cor- 
relators decay quickly in the confinement phase. In the 
deconfinement phase the decay is significantly slower for 
all Gu (t). When we switch on a magnetic field with the 
strength qB = (0.63 GeV) , in the confinement phase the 
correlator G zz (r) decays much slower and is significantly 
larger than zero for all r, much like in the deconfinement 
phase. In contrast, the correlators for the perpendicu- 
lar components of the current G xx (r) and G yy (r) decay 
somewhat quicker than in the zero field case. In the 
deconfinement phase all the correlators are practically 
unaffected by the magnetic field. 

We now apply the Maximal Entropy Method [10, 15, 
16] to extract the spectral functions (2) from the correla- 
tors (1). Our analysis is similar to that of Refs. [10, 16]. 
We used the model with the default guess m (w) = 
fho(b + aw) [10]. Some spectral functions at different 
temperatures and magnetic fields are plotted on Fig. 2. 

In the confinement phase and in the absence of mag- 
netic field, the spectral function has a distinct peak near 
w w 1 GeV, which corresponds to the mass of the p me- 
son in quenched SU (2) lattice gauge theory [16, 17]. The 
width of this peak in quenched approximation is a lat- 
tice artifact [16], and should decrease for finer and larger 
lattices. The spectral function in the limit of zero fre- 
quency, p^ (0), is equal to zero within error range. This 
indicates that in the absence of an external magnetic field 
the vacuum of quenched QCD is an insulator, in agree- 
ment with the results of [15, 16]. When the external 
magnetic field is applied, the peak grows and the spectral 
function becomes nonzero in the limit of zero frequency. 
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FIG. 1. The correlator (1) in the confinement (left) and in the deconfinement phases (right) at T — 350 MeV. 
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FIG. 2. Spectral functions pij (w) in the confinement and FIG. 3. Electric conductivity of quenched QCD as a function 



deconfinement phases. 



of an external magnetic field at different temperatures. The 
points for <r zz and a xx at T>T C coincide within the errors. 



For other components of pij (w) nothing changes quali- 
tatively, but the peak which corresponds to the p-meson 
becomes somewhat smaller and shifts slightly to larger 
w. The conductivity stays equal to zero within the error 
range. Thus when the external magnetic field is applied 
to the quenched vacuum of the SU (2) lattice gauge the- 
ory, the vacuum acquires nonzero conductivity, but only 
in the direction of the magnetic field. 

In the deconfinement phase at zero magnetic field, the 
spectral function is nonzero at w = and has a smooth 
peak near w ~ 2 GeV. Thus quenched SU (2) lattice 
gauge theory is a conductor above the deconfinement 
phase transition [10, 15]. Since the shape of the cor- 
relator dj (r) is practically unaffected by the magnetic 
field, the spectral function (w) and the conductivity 
<Jij do not depend on the magnetic field. 

The electric conductivity as a function of external 
magnetic field is plotted on Fig. 3 for the confinement 
and deconfinement phases. In the deconfinement phase 
the temperature is T = 350 MeV. The value of the con- 
ductivity was extracted from the value of the spectral 
function at w = using (4). In the confinement phase 



and at zero magnetic field the conductivity is zero within 
the error range. As the magnetic field is turned on, the 
conductivity a zz in the direction of the magnetic field 
grows, while all other components of remain equal to 
zero within error range. In the deconfinement phase the 
conductivity is isotropic and is practically independent 
of the magnetic field. One can not exclude, of course, 
that there is a weak anisotropy, which cannot be seen 
at the small number of configurations that we have. It 
should be also noted that in our simulations the value of 
conductivity a = 15 ± 2 MeV at T = 350 MeV > T c is 
still much smaller than the results obtained in [10, 15] in 
quenched SU (3) lattice gauge theory with light staggered 
fermions. This difference is likely to be an artifact of a 
quenched theory, since in this case different probes of the 
confinement-deconfincmcnt phase transition might give 
different transition temperatures. In particular, while in 
quenched SU (2) lattice gauge theory the Polyakov loop 
goes to zero at T c = 313.(3) MeV [14], the chiral conden- 
sate is not zero above this temperature [19]. The situa- 
tion might be similar for the insulator-conductor transi- 
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tion, which in the quenched case might be replaced by a 
soft crossover with much smaller conductivity at T > T c . 
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FIG. 4. Electric conductivity in the direction of external mag- 
netic field a zz for different lattice parameters. 

The transport coefficients typically have rather strong 
dependence on lattice parameters. To ensure that the 
nonzero conductivity is not a finite-volume artifact, we 
have also performed the simulations at different lattice 
volumes and lattice spacings (see Table I). The values 
of conductivity <j zz for different lattice parameters are 
plotted on Fig. 4. One can see that as we go to finer and 
larger lattices, the conductivity docs not change within 
statistical errors. 

We conclude that a strong magnetic field can induce 
nonzero electric conductivity of the vacuum of quenched 
non-Abelian lattice gauge theory along the direction of 
the field, turning it into an anisotropic conductor. This 
effect may be called "electric rupture facilitated by mag- 
netic field" ; it may originate from the interplay of gluon 
field topology and an increase in the quark zero mode 
density due to the presence of magnetic field. It can be 
interesting to investigate whether there is some critical 
value of the magnetic field at which the conductivity be- 
comes nonzero. Transitions of this type are known in 
condensed- matter physics [20]. In contrast, in the decon- 
finement phase the vacuum is an isotropic conductor, and 
the value of the conductivity is practically independent of 
the magnetic field. Thus, if a strong magnetic field gen- 
erates an electric current via the chiral magnetic effect in 
a CP-odd background, then the sufficiently strong field 
would guarantee that the charge will propagate through 
the media due to finite electric conductivity in the both 
phases. 

Finally, let us comment on possible experimental con- 
sequences of the phenomenon described above. The ex- 
pectation value (jk (x)ji (y)) is related to the polariza- 
tion of soft photons and to the angular distribution of 
soft photons and dilepton pairs emitted in the collision 
process [15, 21]. One can therefore expect an enhance- 
ment of production rate in the direction perpendicular to 
the reaction plane, possibly resulting in a negative elliptic 



flow for soft photons and dilepton pairs. 
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